2.18 Adaptive multidimensional modeling with applications in scientific 
computing 


► network problems (energy distribution, telecom networks, 
queueing problems, . . . ) 

► viston/graphics (shape reconstruction, reflectance distribution, 
video signal filtering, . . . ) 

► metamodelling (microwave devices, material design, 
computational finance, . . . ) 


S. Becuwe, X. Granier, K, in r t Hout, W.-s. lee 3 R. Lehmensiek, 
R.B. Lenin, M. Lukach, R. Pacanowski, P. Poulin and C Schlick 
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Adaptive multidimensional modeling 
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such that r(xj , . . 
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Adaptive multidimensional modeling 


Tn.mfrl, ■ ■ ■ ,Xd) 


Pn,m ( X 1 p ■ ■ ■ < ^d ) 

tin, m ( x 1p ■ ■ ■ i ) 


m 






+f( £) n (e) >0 

P n , m \ f > M TL r m s? O 


(\ a,m( x f 1 x 1 ) >0 [ Pn,m f< ‘ Qn.,m ^ 0 


cWJD 


strictly convex QP 

nonempty interior 



b- Benchmark 

► Mational Institute of Standards and Technology ( ISilST) 
reference dataset 

► 151 observations 
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best £2-3pproximation 

^ 2 . 2 ^^’ ■ ■ ■ ■ x d ' ) — f' Cl ) minimal 

f =o 

* _ 1.6745 - 0.13927x + 0.00260x 2 
l2 ’ 2 “ 1 - 0.00172x + 0.00002x 2 



residuals, <j = 0.16355 



f ^ - f ^ = 2 [ 3c) = 0.9813 => n = 2, m = 2 


1.56271 - 0.13713x + 0.00261x 2 
1 - 0.00173x + 0.00002x 2 
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f£ J - f<° = 2(2 a) = 0.6542 =* n = 2, m = 2 


1.16217 - 0,1080x + 0.00224x 2 
1 - 0.00223x + 0.00002x 2 




compare to best t^-approximation 


max 

£=G,...,s 


2,2 


( x i 


(£) 


At) 


) — f (£ ) 


minimal 


1.15538 - 0.10751x + 0.00223x 2 
1 - 0.00223x + 0.00002x 2 



residuals, max = 0.3244 
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compare to best too-approximation 


max 

i= o,_ ( s 


r OO (yU) 
r 2,2* x l i 


-f" 1 


minimal 


1.15538 - 0.10751x + 0.00223x 2 
1 - 0.00223x + 0.00002x 2 


os 

rN 



■ 0 * 

-OS 


0 So im ISU im 3 M JDU 3X 400 

residuals 



f>^ - f<' = 2( 1.75a) = 0.5724 =£• r 2 2 and r^' 2 do not satisfy 

=> n = 3, m = 2 
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& Grid structured data 


Ideal iowpass filter: 


H (e lL \ e lt2 ) = 

(ti,t 2 ) € P C [— 7t, 7t] X [—71,71], 

(ti,t 2 ) i P. 


In practice: 

► passband [1 — 6i, 1 + 61], (ti,t2]cP 

► stopband [-62,62], (ti,t 2 ) £ PUT 

► transition band [—62,1 + 61], [ 1 1 , t 2 ) €T 




Grid structured data 


Examples of passband: s + 1 = 33 x 33 




centra symmetric filter 


fen filter 
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& Grid structured data 


Rational models for the parameters 8 % = 0.01 and 82 = 0.02 



1*19.20 Uir *2 3 



1*12,14 (tlr 't2 ) 




Point valued data 


parameters — > physical model — > behaviour 
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& Point valued data 


parameters — > physical model — > behaviour 


| simplify j 




Point valued data 


parameters — > physical model — > behaviour 


| simplify j 


parameters — > metamodel — + behaviour 
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b- Point valued data 


Model of the stripline characteristic impedance Zo(q, e r ) 



Z 0 (q,e T ) 



ri4,io[q,£r) 


Support Pcmts lor P] 

2Sft < - — 

ffl ' + 

* 

15 - * 

10 • 

n a * 

iw ^ .-II 

0.2 04 06 0.8 

s I 1 = 25 data points 



b- Point valued data 

Model of the transmission coefficient S2i(f T vv) of two inductive 
posts in rectangular waveguide 





argfSaiff.w)) arg(r m o(fr w)) 


re 

r$ 

I Ia 
f 10 
a 
e 

T 6 * Id ti ti 

rroifal 

$ + 1 =40 data points 
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Interval valued data 


A European call option gives its holder the right (but 
not the obligation) to purchase from the writer a pre- 
scribed asset for a prescribed price at a prescribed time 
in the future. 


T expiry date (0 ^ t < T) 

E strike or exercise price 

S asset price S t ^ 0 

r annual interest rate (constant) 

a market volatility 





Interval valued data 


Black-Scholes PDE 


8C 1 2c2 9 2 c dC 

k + 2 gS asr + TS as— rC = 0 


C(S, T) = max(S — E, 0) 
C(0,t) = 0, O^tsST 

C(S, t) fa S, large S 
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Interval valued data 


Typically a few million values computed: 

(s (w) ,t tl ' k) ) , (), k) = O(10 4 ) 

* fit model through subset of s + 1 data points 

► check model on all available values 

► increase s and update model 




Interval valued data 


Graph in (u, v): u. = InfS) — Inf H) 4- rt, v = cr \/t 

C Sl‘n,m ( U, v) 


n — 16, 


nr = 20, f^ £) -f<’ = 0.005 
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Interval valued data 




ri 6 , 2 o[u,v) 


relative error 


, m ( "Li* V ) 


degree 4 + aisti 5 + a^v 5 
degree 5 




Scattered interval data 


The Bidirectional Reflectance Distribution Function 
p(0t T tt>i,0v,4^v) describes how a materia! reflects light 
from surfaces. 


1 lighting direction 

v viewing direction 

0 zenithal angle 

cj> azimuthal angle 



chrome steel 


fabric beige 




325 



Scattered interval data 


For isotropic materials 90 x 90 x 180 w 1.45 million measured 
BRDF samples (RGB values): 

► p(0i t 4 h, 0 v , 4>v) ft* t n m (0i u 0^ ) 

► a priori error control (3-5%) on all data points 1*12 Mb) 

► a posteriori error control on all measured samples 




Scattered interval data 


Rendered example: blue-metallic paint 



Original (33MB) Approximation (1.15KB) 
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I* Sparsity 


Compressive sensing recovers a K-sparse signal from only 
M ^ K measurements without loss of information. 


x(t) = 2cos(5t) — 15cos(14t] + cos(26t) 0 ^ t ^ 2 tt 
4- 5cos(35t) + noise ([— 0.1, 0.1]), 





Sparsity 


tj = j = o, . . . , 7, s + 1 = 8 
Choice of data points allows to recover the frequencies first, 

5 14 26 35 

and afterwards the coefficients by fitting the same data, 
2.004 - 14.93 0.9668 5.007 

x(t) = 2cos(5t) - 15cos(14t) +cos(26t) +5cos(35t) 
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